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ABSTRACT 

We report the discovery of the first high-amplitude 8 Scuti star in an eclipsing 
binary, which we have designated UNSW-V-500. The system is an Algol-type 
semi-detached eclipsing binary of maximum brightness V = 12.52 mag. A best- 
fitting solution to the binary light curve and two radial velocity curves is derived 
using the Wilson-Devinney code. We identify a late A spectral type primary 
component of mass 1.49 ±0.02 M Q and a late K spectral type secondary of mass 
0.33 ± 0.02 M , with an inclination of 86.5 ± 1.0°, and a period of 5.3504751 ± 
0.0000006 d. A Fourier analysis of the residuals from this solution is performed 
using PERIOD04 to investigate the 5 Scuti pulsations. We detect a single pulsation 
frequency of fi — 13.621 ± 0.015 cd _1 , and it appears this is the first overtone 
radial mode frequency. This system provides the first opportunity to measure 
the dynamical mass for a star of this variable type; previously masses have been 
derived from stellar evolution and pulsation models. 
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1 INTRODUCTION 

The serendipitous arrangement of a pulsating star in an 
eclipsing binary system represents a unique laboratory 
for astrophysical measurements. The binarity constrains 
the physical and geometrical parameters of the system, 
and can also assist in mode identification in the pulsa- 
tions. Since the pulsating star and non-pulsating com- 
panion can reasonably be assumed to have formed from 
the same parent cloud, we can utilise information from 
the non-pulsating companion in identifying stellar evolu- 
tion models for pulsating stars. We present here the first 
known example of a high-amplitude 8 Scuti (HADS) star 
in an eclipsing binary system, designated UNSW-V-500. 

5 Scuti stars are the main-sequence analogues of 
Cepheid variables. They are late A to early F spec- 
tral type, and pulsate with periods of between 1-6 
hours. They are typically on or slightly above the main- 
sequence. Low-amplitude 8 Scuti stars typically pulsate 
in many higher, non-radial modes simultaneously with 
amplitudes of less than 0.05 mag. High-amplitude 8 
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Scuti pulsate primarily in the radial modes and have 
higher amplitudes, with the conventional cut-off given as 
A\r ^ 0.30 mag. Fig.[T]shows the 8 Scuti region of the HR 
diagram; high-amplitude 5 Scuti stars are constrained to 
a narrower range in T c ff of width ~300 K within this re- 
gion. The subset of low metallicity Population II HADS 
stars are designated as SX Phe stars. Prior to the dis- 
covery of UNSW-V-500, all bright field HADS stars were 
identified as pulsating in the fundamental radial mode. A 
significant fraction (~40 per cent) are double-mode pul- 
sators, additionally pulsa ting in the radial first overtone 
mode (McNamara 2000a). However, SX Phe stars have 
been identified in globular clusters as pulsatin g in the first 
overtone rather than the fundamental mode (|McNamaral 
l2000bl : lNemec et al.|[l994h . 



ISovdugan et alj l|2006l ) present a list of 25 confirmed 
eclipsing binary systems with pulsating components in 
the 8 Scuti region of the instability strip. All of these 
have low pulsation amplitudes, ranging from A\r = 0.007- 
0.02 mag and up to Ab = 0.06 mag. Thr ee systems 
that have been s tudied in detail a r e Y C am (|Kim et al.l 
I2002D. AS Eri ilMkrtichian et al.l I2004D and AB Cas 
(|Rodri'guez et allboCMl h" Parameters of these systems are 
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Figure 1. A HR diagram showing the positions of 5 Scuti 
stars in ec lipsing binaries. The s olid triangles are data from 
Table 4 of ISovdugan et alj (2006). The open circles are high- 
amplitude S Scuti stars from Table 2 of lMcNamaral (|2000bh . 
UNSW-V-500 is shown as an open square. The observational 
red and blue edges of the classical instability strip are shown, 
as well as the theoretical zero-age main sequence. 



shown in Table [T] with the parameters of UNSW-V-500 
shown for comparison. 

We note that although 35 8 Scuti sta rs in eclipsing bi- 
nary s ystems have be en identified (25 in Sovduga n et alj 
(|200fJ ); 9 recently in IPigulski fc Michalskal <|2007f ) : and 
the subject of this paper), UNSW-V-500 is the first 
HADS star. This is a curious statistic, since with larger 
amplitudes, surveys might be expected to be observation- 
ally biased towards finding HADS stars in these systems. 
A census of the Rodriguez, Lopex-Gonzalez & Lopez de 
Coca (2000) catalogue shows that the detected fraction 
of HADS stars is ~25 per cent of the total S Scuti popula- 
tion. One assumes that this fraction is influenced by the 
same observational selection bias affecting the detection 
of HADS stars in eclipsing binary systems; we are seeing 
a detection rate nearly an order of magnitude below this. 
However, it is difficult to draw conclusions from the small 
numbers of binary systems that have been found. 



2 OBSERVATIONS 
2.1 Photometry 

UNSW-V-500 was initially observed over 29 nights from 
February 2006 to April 2006. The photometric /-band 
observations were performed with the 0.5-m Automated 
Patrol Telescope at Siding Spring Observatory, Australia. 
The observations formed part of an extrasolar planet 
transit search being undertaken b y the University of New 
South Wales l|Hidas et al.ll2005l ). The CCD has 770 x 
1150 pixels and images a 2 x 3 deg 2 field with a relatively 



low spatial resolution of 9.4 arcsec per pixel. We have 
used a customised aperture photometry data reduction 
pipeline to construct our target light curves. The results 
of the project include the disco very of a new eclip sing 
system of K7 dwarf components l|Young et al.|[2006l ). 

Identification of UNSW-V-500 as demonstrating 
both eclipsing and pulsating variations was made during 
routine cataloguing of the variable light curves detected 
in the transit search, to be published separately. The pho- 
tometry aperture used in the reduction process, which is 
nearly 1 arcmin in diameter, can usually be expected to 
contain more than one star due to crowding effects in the 
target fields (typically chosen to have a Galactic latitude 
of 10-20°). Higher resolution images from the Digitized 
Sky Survey catalogue showed that in the case of UNSW- 
V-500 the photometry aperture contained one bright cen- 
tral star and 6 additional stars at least 3.5 magnitudes 
fainter. Due to the large amplitude of the S Scuti pul- 
sations (diluted to ~0.1 mag in the original photometry 
aperture) the system was identified with the bright cen- 
tral star, aj2000 = 13h 10 ' 18 " 7 > 5 J2000 = ~ 45 ° 9 ' 13 "- A 
catalogue search revealed that this system had been pre- 
viously observed and identified in the A ll Sky Automated 
Survey Catalog of Variable Stars III (|Poimanskil l2004h 
as an eclipsing binary system, designated ASAS131018- 
4509.2. From their data they measured an initial epoch 
of T = JD2451892.6 and a period of P = 5.350479 d. 
However, their precision was insufficient for detection of 
the 5 Scuti pulsation. 

In the original run of 29 nights one secondary min- 
imum and two partial secondary eclipses were observed, 
but only two primary eclipse egresses. In order to improve 
coverage of this part of the light curve, UNSW-V-500 was 
observed again in the same configuration as described 
previously on two nights in February 2007 and March 
2007 at the predicted times of primary eclipse. These two 
nights and the 22 best nights of the 2006 data are shown 
in Fig. [2] phased at a period of 5.3504751 d. The primary 
eclipse data conclusively confirmed that the light curve 
was not the result of a 5 Scuti star blended with a back- 
ground eclipsing binary: the flat-bottomed eclipses show 
no sign of the pulsations that are evident in the remain- 
der of the light curve. Therefore, the <5 Scuti star must be 
fully eclipsed by the secondary component of the binary 
during primary eclipse. The inset in Fig.[2]shows the flat 
primary eclipse in more detail. 

In order to confirm the identification of the bright 
central star in the original photometry aperture as the 
eclipsing binary, higher spatial resolution observations 
were obtained with the 40-inch telescope at Siding Spring 
Observatory. The WFI CCD mosaic was used, with an 
image scale of 0.38 arcsec per pixel. Observations were 
taken on a single night in January 2007, in the Johnson V 
filter. These data were reduced using a modified version 
of our aperture photometry pipeline. The identification of 
the pulsating star was confirmed and the data are shown 
in Fig. [3] as solid circles. For comparison, the light curve 
of a nearby star of similar magnitude, GSC0824700373, 
is shown as open squares. The lower limit on the ampli- 
tude of the 8 Scuti pulsation, diluted in these data by 
light from the secondary, is Ay = 0.21 ± 0.02 mag. Com- 
bined with a primary eclipse depth in the / band of ~60 
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Table 1. A sample of 8 Scuti stars in eclipsing binary systems. Ay is the amplitude of the 8 Scuti pulsations. 



Name 


V 


A v 


Spectral Type 


Ppuls 


Porb 


Inclination 




mag 


mag 




(d) 


(d) 


o 


Y Cam 


10.56 


0.04 


A7V 


0.063 


3.3055 


86 


AB Cas 


10.16 


0.05 


A3V 


0.058 


1.3669 


87.5 


AS Eri 


8.31 


0.0068 


A3V 


0.016 


2.6642 




UNSW-V-500 


12.52 a 


~0.35±0.5 


A7V 


0.0734±0.0001 


5.3504751±0.0000006 


86.5±1.0 



a V m ax from the ASAS catalogue. 





4110.2 -1110.23 
HJD 2450000.0+ 



Figure 3. Light Curves observed with the higher spatial res- 
olution of the 40" telescope. The solid circles show the 8 
Scuti pulsation of our target, and the open squares show the 
light curve of the nearby star GSC0824700373 for comparison. 



Figure 2. The phased light curve of 24 nights of data taken 
with the Automated Patrol Telescope in the I band. The upper 
curve is the original data, with the scatter outside the primary 
eclipse due to the S Scuti pulsation. The lower curve is the 
same data with the 5 Scuti pulsation reconstructed from the 
frequency analysis and removed. In both cases, the solid line is 
the fit to the original curve using the Wilson-Devinney code. 
Panel (a) shows the primary eclipse in more detail — there is 
no evidence of <5 Scuti variations in this region. 



per cent, the final lower limit is Av = 0.35 ± 0.05 mag, 
confirming this star as a high-amplitude 8 Scuti star. 

2.2 Spectroscopy 

Several medium resolution spectra (R ~ 6000) were ob- 
tained over two nights in February 2007 with the Double- 
Beam Spectrograph on the 2.3-m telescope at Siding 
Spring Observatory. The wavelength range covered was 
3900-4400 A in the blue arm and 8000-8900 A in the red. 
The spectra were reduced using standard IRApQ spec- 

1 IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universi- 
ties for Research in Astronomy, Inc., under cooperative agree- 
ment with the National Science Foundation. 



troscopy routines. The observations were alternated with 
arc spectra of FeAr in the blue and NeAr in the red. The 
flux calibration of the system was performed using the 
standard stars HR4469 and HR4963. The spectra were 
rebinned to a resolution of 10 A using the IRAF rou- 
tine REBIN and compared with UVILIB template spectra 
(|Pickleslll998l ). Visual inspection resulted in the classifi- 
cation of the spectra as an A7V star. The phase coverage 
was insufficient to measure the dynamical mass. 

To obtain sufficient phase coverage, we obtained 
additional spectra with the same instrument on five 
nights in May 2007. The gratings and grating angle were 
changed to give a wavelength coverage of 3600-4700 A m 
the blue arm and 6000-7000 A in the red. The same 
procedure of alternating observations with arc spectra 
for wavelength calibration was followed, however to im- 
prove the stability of the calibration between arc spec- 
tra we used the night sky lines present in each spec- 
trum in the red half of the data for additional calibra- 
tion. The data were then continuum normalised. The red 
data clearly show spectral features of both the primary 
and secondary components, and also a significant compo- 
nent of Ha emission. Given the semi-detached nature of 
UNSW-V-500 (see Sec. [3]) this may indicate the existence 
of a gas stream between the two components (see for in- 
stance the_set_of_well_ol2se^ 

ries in lRichards fc AlbrightJ (|l999l ); IVesper et all (|200lh V 
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The blue data are entirely dominated by the spectrum 
of the primary component. Therefore the blue data were 
used to identify the spectral type of the primary, by using 
the prelimary identification from the earlie r data and the 
ATL AS9 synthetic stellar template library (jMunari et al] 
2005), rebinned as previously. We performed a least- 
squares fit and identified the T e g = 7500.R', log g — 4.0, 
[Fe/H]= —0.5 template as the best fit. From the residuals 
to this fit, we attempted to match the spectrum of the 
secondary component. A preliminary light curve analysis 
had indicated a secondary component with a temparaturc 
of ~4200K; hence the least-squares fit to the residuals was 
restricted to the ATLAS9 models with T c « = 4250-fY and 
[Fe/H]= —0.5, since we can assume the binary system 
will have a common origin and thus metallicity. The best 
fit was achieved with the log g — 3.0 template. 

We note again that the pulsating primary compo- 
nent is essentially fully eclipsed by the secondary star at 
primary eclipse. Therefore, a high resolution spectrum 
during the time of primary eclipse would necessarily be a 
spectrum of the secondary star, and would be useful for 
constraining the stellar spectral type and physical param- 
eters derived via light curve-fitting in Section [3] 
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Figure 4. The radial velocities of the primary (solid triangles) 
and secondary (open circles) components. The lines are the 
best fitting sine curves. 



2.3 Radial velocity analysis 

In order to extract the radial velocities of the 
two b inary components, we used the program TOD- 
COR (jZucker fc Mazehl ll994T ). which performs a two- 
dimensional correlation between two supplied template 
spectra and an object spectrum of a binary system. Using 
the two ATLAS9 templates identified previously, radial 
velocities were extracted for the majority of the spec- 
tra we had obtained. The correlation was limited to the 
wavelength region 6200-6530 A to avoid the Ha emis- 
sion noted previously. The flux ratio of the secondary to 
the primary template spectra was left as a free parameter 
for TODCOR, and was found to vary from 0.2-0.4 with 
phase. 

The radial velocities are shown in Fig. [4] with the 
primary component shown as circles and the secondary 
component as squares. The lines are the best-fitting sine 
curves, with velocity amplitudes of K\ = 27.0 ± 1.8 km 
s _1 and K2 = 121. 2± 1.4 km s" 1 , indicating a mass ratio 
of 0.22±0.02, and a systemic velocity of 43.6 ± 0.9 km 
s" 1 . 

The large scatter in the primary component data of 
~20 km s _1 is due to the S Scuti radial velocity pulsa- 
tions, and is similar to the radial velocity amplitude of 
other HA PS stars we have m easured with the same in- 
strument l|Derekas et al] 12006). The frequency spectrum 
of these data were analysed in the same manner described 
in Sec. |4] and two peaks corresponding to the orbital and 
pulsational periods (5.36 d and 0.073 d) were identified, 
a second confirmation that this is not a blended system. 

3 BINARY SYSTEM 

In order to fit the orbital parameters of UNSW-V-500, we 
applied the Wilson-Dev inney code l|Wilson fc Devinnevl 
ll97ll ; IWilson1ll979l . ll990l ) to the APT light curve and the 



two radial velocity curves simultaneously. This system 
is an Algol-type semi-detached eclipsing binary system, 
with the secondary star filling its Roche lobe, and con- 
sequently the code was operated in mode 5. The effec- 
tive temperature of the primary was fixed at T\ — 7500K 
from the template fit. The gravity brightening coefficients 
were set to 1.00 for the radiative primary component and 
0.32 for the convective secondary component. The albe- 
dos were set to the standard theoretical value of 1.00. 
The bolometric and bandpass-specific limb darkening co- 
ef ficients were adopted from values for the closest models 
in Ivan Hammd (|l993n . The third light Z3 was assumed 
to be non-zero due to the crowded photometry aperture, 
and was allowed to vary as a free parameter. It was given 
an initial value of Z3 = 0.2 from an estimate of the max- 
imum total contribution to the normalised flux at phase 
0.25. The eccentricity was assumed to be ~0 due to the 
secondary eclipse occurring at a phase of 0.5. To confirm 
this, the eccentricity was allowed to vary and did not 
result in any significant improvement in the fit, so was 
fixed at for the subsequent fitting. The semi-major axis 
was fixed at 15.69 Rq from the total mass and period of 
the system, and the systemic velocity was fixed at 43.6 
km s _1 . The mass ratio was fixed at 0.22 from the radial 
velocity data. The free parameters were thus the third 
light £3, the inclination i, the effective tem perature of th e 
secondary T2, the potential (as defined bv lKopal l)l954h ) 
of the primary fii and the luminosity of the primary L\ . 
The values of these parameters used in the final solution 
are shown in Table 2, as are the derived quantities for 
the two components of mass, radius, log g and M(, (. 

The solid line in Fig. [2] shows the best-fitting solu- 
tion. The high inclination (i — 86.5 ± 1.0°) is indicated 
by the flat bottom of the primary eclipse. In fact this 
is the first eclipsing binary system containing a pulsat- 
ing component to demonstrate a flat-bottomed eclipse, 
with the possible exception of the r ecent discovery of the 
pulsating component of HD 99612 (Pig ulski fc Michalskal 
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Figure 5. The residuals in the original data after the subtraction of the best-fitting binary solution, between phases 0.1 and 0.9. 



Table 2. The parameters for the binary system solution. 
Starred quantities indicate the free parameters in the Wilson- 
Devinney code. Quoted errors for the starred quantities and 
the calculated quantities (L 2 and O2) are standard deviations 
produced by the Wilson-Devinney code. a These are the band- 
pass luminosites in the /-band. b This is the corrected value 
of the third light for reference phase 0.25. 



Parameter 


Value 


e 


0.0000 


q* 


0.22±0.02 


Ti (K) 


7500 


T 2 * (K) 


3850±20 


Li* (L Q r 


6.96±0.03 




3.89±0.03 


Qi* 


6.9±0.1 


n 2 


2.3±0.1 


i* (°) 


86.5±1.0 


h* b 


0.096±0.005 


Mi (M s ) 


1.49±0.02 


M 2 (M s ) 


0.33±0.02 


Ri (R@) 


2.35±0.02 


R 2 (R e ) 


4.04±0.01 


M bolil 


1.80±0.02 


M boU2 


3.53±0.02 


log 31 


3.87±0.01 


log 92 


2.74±0.01 



2007). The normalised third light contribution is found 
to be 0.096 at the reference phase of 0.25. 

We note that the log g values that have been derived 
in the Wilson-Devinney fit (3.87 ± 0.01 and 2.74 ± 0.01 
for the primary and secondary components respectively) 
confirm the estimated values from the synthetic template 
fitting (4.0 ± 0.5 and 3.0 ± 0.5) in Sec. H2] 

From the T^ff and derived mass we have identified 
the two components as a late A spectral type primary, 



confirming the A7V classification, and a late K spec- 
tral type secondary. Using the derived parameters, we 
attempted to fit the positions of the two components in 
the HR diagram with the Y 2 isochrones l|Yi et al.l l200ll ; 
iKim et al1 l2002) in order to find an age estimate for the 
system. However, we found that the isochrones and evo- 
lutionary tracks were unable to reproduce the current po- 
sitions of the components; hence we conclude that there 
has been significant mass transfer to the pulsating pri- 
mary component from the secondary component. This 
component appears in a much more evolved state despite 
its lower mass, a well-recognized phenomenon known as 
the Algol paradox. We caution that standard evolution 
and pulsation models for HADS stars may not apply 
to UNSW-V-500 due to its binary evolution. However, 
we do note that it is well described by the ML3 mass- 
luminosity relation for HADS stars sho wn in fig. 1 of 
IPetersen fc Christensen-Dalsgaardl (|l996l ) , which is based 
on models with a metal content of Z = 0.02. 



4 PULSATION 

Once the binary solution has been subtracted, an analysis 
of the S Scuti pulsation can be performed. We have used 
data from phase 0.1 to 0.9 for this analysis, discarding 
those data around the primary eclipse where the HADS 
star is completely eclipsed. This reduces the total number 
of data points by 700, or ~25 per cent. The residuals from 
the binary subtraction between phase 0.1 to 0.9 are shown 
in Fig. 

Fig. [6] shows the freque ncy analysis, as perfo rmed 
with the program PERIOD04 (|Lenz fc Bregerll2005l ). The 
spectral window is shown in panel (a). Panel (b) shows 
the initial periodogram. The dominant frequency is iden- 
tified as /1 = 13.621 ± 0.015 cd _1 , i.e. a period of 
0.0734 ± 0.0001 d, which is typical for 5 Scuti stars. The 
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Figure 6. The frequency analysis of the pulsation. Panel (a) 
shows the spectral window. Panel (b) shows the strongest fre- 
quency /i at 13.621 cd — 1 . Panel (c) shows the frequency spec- 
trum with /i removed. On this scale, /4 is coincident with fo. 



5 SUMMARY 

We have presented here the detection of the first exam- 
ple of a high-amplitude 8 Scuti star in an eclipsing binary 
system, and the probable first detection of a field HADS 
star pulsating in the single first overtone radial mode. 
Several HADS stars have been det ected in binary sys - 
tems previously (including SZ Lyn (Derekas ct al. 2003) 
with a period of 1190 d; and RS Gru Uoner fe Lanevl 
|2004| ) with a period of approximately 2 weeks), however 
these are much wider systems. This new fully eclipsing bi- 
nary opens up many further opportunities for studies of 
HADS stars and pulsating stars in binary systems. Many 
of the poorly understood processes governing the effects 
of mass transfer, tidal interactions, rotation, convection 
and magnetism on 5 Scuti pulsations may be explored 
with this system. 
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four addition al frequencies identi fied with a S/N > 4.0, as 
suggested bv lBreger et al.l (| 19931 ). are shown in panel (c), 
after removal of /i . These can be identified as low-power 
frequencies, probably due to artefacts of the binary sub- 
traction (/ 2 = 0.187 ±0.033 cd" 1 and / 4 = 0.255 ±0.084 
cd" 1 ), or harmonics of /i (/ 3 = 2/i = 27.242 ± 0.084 
cd -1 and f s = 3/i = 40.86 ± 0.14 cd" 1 ). The absence of 
any additional frequencies in the S Scuti range supports 
the identification of a HADS star oscillating in a single 
radial mode to the limits of our detection. 

As an additional check, these frequencies were re- 
moved from the original data. The resulting light curve 
is shown in Fig. [2] offset below the original data. The 
Wilson-Devinney code was re-run using this second light 
curve, with no significant change in the derived parame- 
ters. 

Using the relation from lBreger et all l|l993l ). we can 
calculate the pulsation constant Qobs- 

logQobs = -6.456 ± logP + 0.51og<? + 0.1M bol + logT e (1) 

We find a Q obs of 0.025±0.001. 

iPetersen fe Jorgensenl (1 19721 ) derived theoretical pul- 
sation constants of Qo — 0.0333 for the fundamental 
radial mode, Q\ — 0.0252 for the first overtone radial 
mode, and Q2 = 0.0201 for the second overtone radial 
mode. Inspection of the HADS stars catalogued in 
iMcNamaral (|2000al ) reveals that of the 26 well studied 
field stars, all have observed pulsation constants in the 
range 0.0309-0.331, indicating they are pulsating in 
the fundamental mode. UNSW-V-500 appears to be 
the first identified to be pulsating primarily in the first 
overtone radial mode, joining a number of SX Phe stars 
in globular clusters to have been identified in this mode 
<|McNamarall2000bl) . 
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